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ABSTRACT 


The  results  of  coherence  measurements  of  three  free  precession 
proton  magnetometers  are  discussed.  The  measuring  and  data  process 
sing  instrumentation  is  described.  Ihe  noise  associated  •with  the 
instrumen'bation  has  been  measured  and  presented  in  statistical  fom. 
The  results  of  measurements  with  elementary  local  perturbations  are 
included  to  illustrate  the  target  detection  problm.  The  limitations 
of  the  instrument  for  "target  detection  are  discussed  and  imporvements 
suggested. 

The  ■writer  ’rfishes  to  express  appreciation  to  Professor  Carl  E. 
Menneken  for  suggesting  the  problems,  and  Professor  Harold  A.  Titus 
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1,  Introduction 


The  use  of  magnetometers  for  target  detection  requires  identic 
fication  of  signals  in  the  presence  of  the  noise  associated  with  the 
earth* s  magnetic  field.  The  signal,  in  mangr  cases,  will  be  smaller 
in  magnitude  than  the  normal  diurnal  variation  of  the  field,  local 
fixed  perturbations  and  geomagnetic  piilsations.  Ihe  system  noise 
is  of  major  concern  in  a  detection  scheme  using  magnetometers.  This 
paper  presents  the  preliminary  results  of  a  stucfy  of  three  proton 
free  precession  magnetometers  operated  under  non  laboratory  condi=> 
tions  at  Naval  Postgraduate  School  La  Mesa  Village  site. 

The  free  precession  magnetometer  consists  of  a  coil  of  wire 
in  which  a  source  of  protons  is  placed.  A  strong  field  is  applied 
to  orientate  the  proton  spins  in  a  direction  approximately  nomal 
to  the  earth's  field.  The  polarizing  field  is  then  removed  in  a  man- 
ner  such  that  the  spins  cannot  follow  the  field.  The  earth's  field 
then  causes  the  protons  to  process  about  the  lines  of  magnetic  force 
at  the  Lamour  frequency  given  by 2 


where  '^p  s  gyro  magnetic  ratio 

s  magnetic  raoment/angular  momentum  (1) 


and 

|-j  =  magnitude  of  external  field 

(earth's  field) 

The  net  precessing  magnetic  moment  causes  a  voltage  to  be  induced  in 
a  coil  surrounding  the  sample.  Since  the  processional  frequency  is 
dependent  only  on  the  external  magnetic  field,  the  measurement  of  H 
consists  of  measuring  the  frequency  of  the  induced  emf.  The  inherent 
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accuracy  of  the  free  precession  magnetometer  is  determined  by  the 
spectral  line  width  (l)s  which  is  related  to  the  relaxation  time 
T2>  and  is  given  bys 


where 


L 


=  3.t 

H  s  Half  intensity  width 
T2  s  1  sec  (meas) 
s  10°^  GAUSS 


(2) 


The  gyromagnetic  ratio  has  been  raeasxired  by  the  NBS  to  better  than 
1  part  in  10^,  Full  advantage  of  this  inherent  accuracy  has  not  been 
utilized  because  of  the  problems  associated  with  the  frequency  measure- 
ment.  The  frequency  measurement  is  complicated  by  the  fact  that  only 
a  limited  time  is  available  for  measurement.  The  signal  received 
from  the  magnetometer  is  in  the  few  microvolt  rangej  hence  locally 
generated  electronic  noise  and  transients  from  the  polarizing  opera¬ 
tion  cause  the  signal  to  noise  ratio  to  be  small.  Standard  counting 
techniques  are  used  for  measurement.  To  get  the  required  accuracyj 
the  frequency  must  be  obtained  from  period  or  multiple  period  average 
measurements.  Ihe  process  is  digital  in  natxire  and  lends  itself 
very  nicely  to  computer  processing  of  the  data. 

The  main  disadvantage  is  that  any  reasonable  attempt  to  obtain 
an  analog  signal  proportional  to  the  measurement  is  an  analog  of 
period  rather  than  frequency?  therefore,  the  analog  is  the  inverse 
of  the  quantity  of  interest.  This  does  not  present  serious  problems 
for  target  detection  since  the  absolute  magnitude  of  field  is  not 
of  interest. 
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2.  Ins  trumenta  tion 


2,1  Measurement 

The  engineering  aspect  of  the  measurement  problem  can  be  stated 
simply  as  follows? 

Given?  Three  transient  type  signals  of  frequency  approximately  two 
kc/s  occuring  at  the  same  time 5  with  usable  duration  of  about  3/U 
sec.  The  signals  are  in  the  microvolt  range.  The  recurrance  rate 
is  to  be  2  sec.  The  measuirements  must  be  conducted  in  the  presense 
of  high  electronic  noise  such  as  relay  transients^  contact  noise 
and  transients  from  the  collapsing  polarizing  field. 

Find?  The  frequency  of  the  oscillations  accurate  to  1  part  in  10^ 
or  better,  and  suitable  means  of  data  aquisition  and  processing. 

The  funds  available  excluded  the  possibility  of  using  commer^ 
cially  constructed  instruments  or  components.  Also,  modifications 
and  interfacing  wovild  have  to  be  done  if  components  were  used,  Hsnce, 
the  basic  measuring  instrumentation  was  constructed  at  NPGS,  No 
detailed  circuit  design  is  presented|  however  the  circuit  diagrams 
are  included  in  Appendix  IV, 

Assume  the  precession  signal  has  been  amplified  in  a  suitable, 
low  noise,  narrow  band  preamplifier,  and  converted  to  a  signal  hav= 
ing  fast  waveform.  Referring  to  the  simplified  block  diagram  of 
Fig,  1,  the  basic  counter  operation  is  as  follows?  Let  the  slow 
and  fast  counters  consist  of  flip=flop  counting  units,  connected  to 
form  a  binary  counter.  The  gate  control  flip=flop  is  capable  of 
being  triggered  into  either  state.  The  multiplier  converts  the  clock 
frequency  to  a  harmonic  of  the  precision  100  kc/s  oscillator. 
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BASIC  COUMTER  UNIT 


Input 


100  kc 

Ctysta.1  Oscillator 


Figttre  1 

Let  all  counters  be  reset  and  ready  for  operation.  The  unknown 
input  signal  is  applied  to  the  gate.  After  a  delay  of  approximately 
1^0  msj  an  external  signal  triggers  the  gate  control  to  open  the  sig=* 
nal  gate.  The  long  delay  is  necessary  to  allow  for  damping  the  polari= 
zing  field  and  relay  transients.  The  first  pulse  through  the  signal 
gate  opens  the  clock  gate.  The  fast  counter  then  counts  the  number 
of  clock  pulses  occuring  during  the  time  interval  determined  by  the 
slow  counter.  The  slow  counter  is  increased  one  for  each  cycle  of 
xinknown  input  signal.  After  the  slow  counter  has  reached  a  predeter^ 
mined  number,  both  signal  and  clock  gates  are  closed,  and  the  niamber 
in  the  fast  counter  is  a  value  of  the  tine  interval  of  2®  cycles  of 
unknown  input  signal,  or  can  be  reduced  to  the  average  period  of  the 
signal. 
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Let  us  define 

N  s  number  in  fast  counter 

2*^  s  number  of  cycles  of  unknown  signal 

tp  g  period  of  clock  (^^sec.) 

tjj  g  period  of  unknown  signal 


Then  we  have 


tn-  N  'tf. 


=  NtpjZ"' 

Substituting  in  (1)  gives 

H  = 

Yp 

and  with  n  s  10, 


(3) 


I  I  ^  ^3‘4^400X  102^4" 

n  -  (li) 

The  accuracy  of  such  a  scheme  depends  ong 

a.  The  accuracy  in  determining  the  zero  crossing  of  the  unknown 
signal. 

b.  Propagation  delay  in  the  slow  counter. 

c.  Gating  speed. 

d«  Accuracy  of  the  clock. 

e.  Plus  or  minus  the  las  bit,  due  to  noncoherent  gating  as 
illustrated  in  Fig.  2. 
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EFFECT  OF  NONCOHERENT  GATING 


clock 


N 


8 


Pulses 


I  I  I  M  M 


N 


1 


Sample  Period—-^ 


Figure_2 


The  errors  due  to  determining  the  zero  crossing  are  usually  the  domi^ 
nant  error  in  using  this  techniquei  however  when  averaged  over  102U 
cycles  they  become  small  and  the  largest  error  is  the  plus  or  minus 
last  bit,  as  shown.  Ihe  accuracy  in  determining  the  zero  crossing 
in  the  presence  of  noise  (2)  isg 


V. 


n 


where? 

t^  g  maximum  error  in  trigger  time  in  sec. 

Vn  s  peak  to  peak  noise  in  volts. 

S  s  slope  of  ideal  signal  at  point  of  trigger  levels 
For  a  sinusoidal  inputs 


S 


Assume  a  signal  to  noise  ratio  2  s  20  db,  and 


W  ^  Z  !C  c  /s 


then? 
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For  102ii  periods  8 

O.  3  X /O  ^seo 

The  plus  or  minus  last  bit  error  if  the  clock  is  100  kCg  is  ^10 
sec.,  and  the  other  errors  mentioned  are  negligible  compared  to  the 
10^,,4^sec«  error. 

Consider  the  bit  limitation  in  detail  for  the  earth’s  field 
50?50  ^  ,  wtiich  is  the  approxinate  value  at  Konterey,  California, 
lets 

N  s  OOI3U05I3  ^  1  bit 
1  bit  10^^,^sec. 

them 

/i  H  gii.i  y 

This  becomes  the  sensitivity  of  the  insiruraentc  If  the  clock  is 
increased  to  1  me /sec,  1  bit  s ands 

^  H  e  0.11  'y'  I  however  other  noise  becomes  appreciable  and 
cannot  be  neglected.  The  circuits  were  designed  for  1  m/see  opera 
tion,  but  the  measurements  to  date  vrere  conducted  at  100  kc/s  and 
300  kc/s  clock  rates.  The  added  complication  is  to  construct  3 
meters  that  id-ll  read  identical,  the  input  circtiitry'  and  slow  count 
ers  must  have  identical  characteristics.  This  is  Impcssible 
a  practical  point  of  view.  The  I'es’alt  is  that  Uie  three  counters 
Kith  an  ideal  signal  input  can  vary-^  1  bit  in  comp-arirjg  their  read 
ings. 

A  simplified  block  diagram  and  photo  of  the  nagnetometer  is 
shown  in  Figs,  3,  4  and  5.  It  must  be  remerfbered  t’nal  operation 
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is  to  be  in  remote  locations  where  extremes  of  temperatures  line 
voltage  fluctuationss  poor  groimdss  etc*  will  be  found|  hence  cir^ 
cuits  that  are  sensitive  to  these  variations  must  be  carefully  com® 
pensated* 

For  the  experiments  under  considerations  locating  the  sensing 
coils  approximately  500*  apart  would  be  sufficientj  however  for  fur® 
ther  experiments  a  distance  of  greater  than  500*  might  be  necessary. 
Therefore,  a  provision  to  place  the  preamps  in  the  line  nearer  the 
sensing  coils  for  improved  signal  to  noise  ratio  was  included.  Bias 
for  the  preamps  is  applied  to  the  signal  line,  from  the  console  to 
eliminate  the  need  for  battery  replacement.  The  relays  that  provide 
the  two  to  three  amps  of  polarizing  current,  and  provide  for  dissapa® 
tion  of  the  stored  energy  when  polarizing  current  is  turned  offs  must 
precede  the  preamp.  Since  this  unit  is  placed  near  the  sense  coil, 
it  must  be  magnetically  clean. 

Referring  to  the  block  diagram  the  components  that  are  common 
to  the  three  units  are  the  clock,  timer,  recorders  and  the  power 
supplies. 
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COUNTER  INSTRUlffiNT 


c 


•H 


On 


40 


'O 

c 

cS 


0) 

E 

Cu 

•H 

§• 

W 


<D  tiO 

3  ^ 
o  p 
o  to 


s 

en 


Pu 

s 

M 


I 

o 


O 


Figure 


11 


2.2  Data  Aquisitton  and  Processir^ 

rw  m  ~P~r^i  I  »i~«  nr  m  rir~B  rriw^iMT^  i^iimiiifr 

Figvire  6  shows  a  block  diagram  of  the  data  collection  and  repro¬ 
duction  system. 

The  scheme  used  for  data  handling  is  certainly  not  optimum  in 
tlie  sense  of  tape  utilization  and  ease  of  operation^  but  an  Ampex 
CPIOO  \m.s  available  and  required  no  additional  expendit-ure  for  tape 
recorder.  Also  vdth  minor  modifications  a  data,  recorder  can  be  uti¬ 
lized  at  a  later  date  if  funds  become  available. 

Considei'  the  form  of  the  data.  The  sensing  coils  are  polarized 
for  one  second,  at  the  anl  of  this  time  the  polarizing  current  is 
turned  off,  the  transients  are  dampedi  and  the  signal  is  fed  to  tlie 
countii^  units.  After  102U  cycles  of  the  input  signal  (about  1/2 
sec.)  the  measurement  is  complete.  In  the  fast  counters,  there  are 
three  binary  nxmibers  to  be  processed.  The  magnitude  of  the  binary 
numbers,  for  H  s  50950  vd.th  clock  frequency  s  100  kc/s,  is  about 
0013l4lU6g  or  l6  bits.  With  the  clock  incx'eased  to  300  kc/s,  the 
number  is  00U25260g  or  18  bits.  Hence,  at  the  end  of  the  counting 
interval,  there  are  1*8  bits  to  be  recorded  in  a  format  foi'  computer- 
input.  The  Avtpex  CPIOO  is  a  "portable"  instrumentation  recorder 
^d.th  lU  parallel  inputs  and  capable  of  speeds  from  1  7/8  ips  to  60 
ips.  All  of  the  record  and  reproduce  amplifiers  that  were  on  hand 
wei*e  not  capable  of  recor'ding  D.C,  levelsj  therefore,  the  information 
was  put  on  tape  in  the  form  of  pulses.  Also,  for  reasonable  spsod 
in  reproduction,  the  data  was  recorded  at  1  ?/6  ips  and  repr-oduced 
at  60  ips.  A  record  of  approximately  6,5  hours  can  be  reproduced  in 
10  minutes.  The  machine  used  to  read  the  raw  data  is  a  CDCI60  com- 
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DATA  AQUISITION  AND  PROCES3IN3 


I 


13 


puter  capable  of  reading  words  of  up  to  12  parallel  bits.  The  mag¬ 
netometer  words  are  broken  into  two,  12  bit  words  and  time  multiplexed 
on  tape  to  satisfy  the  computer  input  conditions.  Hence,  for  each 
sample,  there  are  six  12  bit  words  to  be  processed.  Even  with  the 
time  speed  up  in  reproduction  the  CDC  160  is  operating  below  its 
capabilities.  In  an  effort  to  utilize  the  l60  time,  the  data  was 
put  on  tape,  so  that  when  it  occured,  the  l60  would  read  at  its  maxi¬ 
mum  rate  and  use  the  “no  data  time'*  to  continue  with  preliiainary 
processing.  Channels  13  and  lU  on  the  tape  are  input  ready  and  dis¬ 
connect  signals  used  to  control  the  l60.  And  of  course,  for  reason¬ 
able  machine  time  efficiency,  there  is  on  hand  enough  tape  to  allow 
continuous  recording  for  one  week. 

The  Ampex  CPIOO  was  not  designed  specifically  for  this  type  of 
operation!  hence  the  need  for  an  interface  to  mate  it  with  the  160. 

The  interface  must  provide  the  following  functions g 

a.  Convert  the  pulse  type  data  to  D,C.  level  signals. 

b.  Put  the  data  on  line  when  a  request  is  supplied  from  the 
computer. 

c.  Effectively  disconnect  itself  from  the  lines  when  other 
equipment  is  called. 

d.  Provide  control  information,  i.e.  input  ready  and  discon¬ 
nect  for  the  l60. 

e.  Effectively  isolate  the  computer  from  spurious  transients 
from  the  tape  machine.  These  normally  occur  dui'ing  start  -  stop 
operations,  but  the  recorder  is  somewhat  susceptible  to  power  line 
transients  and  they  sometimes  appear  in  its  output. 
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Figures  8,  9  and  10  shows  the  block  diagram  of  the  interface 
provided*  Fig*  11  is  a  photo  of  the  preliminary  processirig  setup. 
The  cards  were  assembled  and  mounted  as  a  component  on  a  laboratory 
cart,  the  tape  unit  is  placed  on  the  top  of  the  cart  and  computer 
cable  connectors  provided.  This  eliminates  the  need  for  a  permanent 
rack  in  the  computer  laboratoiy. 

The  preliminary  processing  consists  of  the  followings 

a,  Reading  the  raw  data, 

b.  Checking  the  data  and  keeping  a  record  of  the  errors* 
c*  Writing  the  storage  tape* 

d*  Providing  identifying  blocks  on  the  storage  tape. 

e.  Searching  the  storage  tape  for  ary  particular  record  or 
the  end  of  the  data  so  that  one  storage  tape  may  be  used  for  many 
week’s  records, 

f.  Direct  display  of  the  raw  data* 

The  data  storage  format  used  is  illustrated  in  Fig*  7* 


STORAGE  TAPE  FORMAT 
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DATA 
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Figure  7 

With  the  time  base  used  there  are  nearly  380li0  twelve  bit  words 
per  6,5  hours  record.  The  6*5  hours  record  is  one  3600  roll  of  in¬ 
strumentation  tape  at  1  7/8  ips.  Since  the  operation  is  stopped  to 
change  tape,  it  is  convenient  to  put  the  identifying  block  at  the 
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end  of  this  irecord.  Also,  this  number  of  data  words  can  be  process 
sed  in  the  order  of  minutes  by  the  CDCl60li  so  that  if  the  computer 
is  needed  for  other  jobs|  it  can  be  stopped  at  the  identifying  blocks 
and  return  later  to  the  same  point  with  the  proper  index  setting* 

The  block  length  was  set  somewhat  arbitrarily*  The  memory  storage 
limit  of  the  l60,  write  time  of  the  163  tape  unit  and  suitable  re= 
cords  for  handling  in  the  1601;,  are  the  factors  to  be  considered  when 
chosing  the  block  length.  The  instrumentation  recorder  cannot  con-= 
veniently  be  stopped  when  writing  the  storage  tapej  hence  one  data 
point  is  lost  when  transferring  data*  This  represents  22  points 
per  6.5  hours  record  or  less  than  a  minute  in  real  time.  The  block 
length  was  chosen  to  be  3li56]_Q  twelve  bit  xrords  per  block,  the  fig=- 
ure  being  corapatable  xdth  the  above  conditions. 

The  identifying  record  consists  of  xrriting  the  month,  day  and 
year  in  flex  code.  The  remaining  cells  of  the  block  are  loaded  vrith 
zeros.  It  is  then  easy  to  separate  the  data  from  the  identifying 
blocks  since  none  of  the  data  words  will  be  all  zeros. 

The  fact  t’lat  the  data  had  to  be  reconstructed  in  tl  is  Ic'Oi,  x-jas 
also  considered  in  chosing  the  block  length.  The  information  is 
buffered  into  the  l60l;  in  U8  bit  words.  The  blocks  consist  of  an 
even  number  of  data  points,  i.e,  six  12  bit  or  two  li8  bit  words* 

After  a  block  is  read  into  the  machine  and  loaded  into  the  A  register, 
the  magnetometer  readings  can  be  stored  sequentially  in  memory  with 
a  loop  of  a  few  machine  language  instructions.  The  data  analysis 
is  then  done  using  Fortran,  alloxri.ng  easy  and  versatile  manipxalation 
of  the  data.  The  I60  and  I60U  input  programs  are  included  in  Appcn=* 
dices  II  and  III. 
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Reassembly  is  illustrated  below 


160U  WORDS 


Mag  1  (2U  bits) 

Mag  2  (2U  bits)  j 

Mag  3  (2U  bits) 

Mag  1  (2ii  bits) 

Mag  2  (2U  bits) 

Mag  3  (2U  bits) 

Briefly,  the  operation  is  as  follows?  Magnetometer  readings 
are  recorded  for  long  periods  of  time  at  a  remote  site  on  a  protable 
magnetic  tape  units  The  tape  unit  is  then  returned  to  the  process^ 
ing  site,  and  using  a  time  speed  up  of  X  32,  the  intermediate  pro¬ 
cessing  is  performed  and  the  data  stored* 

The  storage  tape  along  with  suitable  programs  is  then  placed 
on  file  for  processing  by  the  computer  facility  at  their  convenience. 
The  tape  unit  is  then  returned  to  magnetometer  site  for  more  records. 
The  new  data  obtained  is  added  to  the  storage  tape  so  as  not  to  de¬ 
stroy  the  existing  data. 
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3.  Data  Analysis 


The  experimental  measurement  program  was  set  up  to  study  the 
coherence  between  the  magnetometersj  i.e,  the  comparison  of  simulta¬ 
neous  samples  of  the  field,  with  the  sensors  located  at  different 
spatial  locations.  The  coherence  was  to  be  examined  with  various 
local  perturbations  in  the  field.  With  the  target  detection  objective 
in  mind,  the  least  complicated  method  of  obtaining  the  target  signal 
will  result  in  the  best  instrumentation  for  field  use.  The  problem 
is  somewhat  analogous  to  the  radar  problem  where  decisions  must  be 
made  in  the  presence  of  noise,  with  associated  probabilities  of  error 
in  detection. 

There  are  several  types  of  geomagnetic  noise  that  occur.  The 
most  dominant  is  the  diurnal  variations.  Superimposed  on  this  are 
smaller  variations  known  as  micropulsations  and  noise  caused  by  lo¬ 
cal  disturbances.  Sudden  eruptions  on  the  sun  sometimes  cause  large 
fluctuations  in  the  field.  Several  studies  of  the  various  noises 
have  been  made  (3),  (Ij.)  and  (^),  For  the  proton  precession  magneto¬ 
meter,  the  instrument  produces  noise  of  major  contribution  to  the 
spectra  of  small  geomagnetic  pulsations.  For  target  detection,  the 
signal  of  interest  will  most  likely  be  of  the  same  order  of  magni¬ 
tude  as  the  small  micropulsations  and  local  distrubances.  For  geo- 
physical  prospecting  and  mapping,  the  signals  will  be  much  larger. 

As  with  any  raw  data,  there  will  be  some  points  that  are  obvi¬ 
ously  in  error.  The  cause  of  errors  is  ustially  uncertain,  but  some 
provision  must  be  made  in  the  processing  to  eliminate  these  errors 
so  as  not  to  distort  the  data.  Care  must  be  taken  when  performing 
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these  operations  so  the  target  signals  are  not  masked.  Records  of 
the  errors  should  be  kept  in  order  that  decisions  about  the  useful- 
ness  of  the  set  of  data  can  be  made. 

The  problem  of  coherence  between  the  spatial  samples  requires 
that  the  instnament  noise  be  determined^  with  an  ideal  signal^  then 
measurements  be  conducted  with  no  local  perturbations  in  the  field, 
and  the  samples  compared.  The  simplest  comparison  is  to  plot  the 
difference  signal  vs.  time.  With  no  local  perturbation,  this  means 
essentially  a  plot  of  the  system  noise  vs.  time|  hence  the  need  for 
filtering  is  obvious.  To  use  this  method  for  target  detection  will 
require  serious  consideration  of  the  filter  time  constants. 

The  types  of  filtering 

1.  A  clipping  filter  (on  the  raw  data)  that  clipped  at  ^  200  y. 
The  number  of  points  outside  of  this  range  are  recorded.  This  range 
was  about  twice  the  expected  diurnal  variation.  The  filter  removes 
the  errors  from  lost  high  order  bits  and  other  instrument  soxu’ces. 

The  typical  error  rat-e  was  about  1/50000. 

2.  A  simple  three  point  average  for  display  of  the  raw  data  on 
the  DD  65  display  xinit  was  used.  This  gives  an  operator  a  fast  vis¬ 
ual  display  of  the  raw  data.  Decisions  can  then  be  made  about  the 
quality  of  the  data. 

3.  Eighteen  point  averaging,  i.e.  36  seconds  real  time,  gave 
partially  smoothed  6,5  hoxirs  graph  records, 

U.  A  simple  single  section  low  pass  filter  was  used  of  the 
form  l/(s  -  a).  The  time  constant  was  set  at  100  sec,  and  20  min, 
in  real  time.  The  implementation  was  done  using  2  transform  (6). 
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The  filter  design  was  done  by  Professor  Harold  A*  Titus 


We  have 


x„(s) 

Itln(S) 


a 


s  -  a 


which  has  a  Z  transform 

X(^(.2)  /  -  e  d 

In  block  diagram  form  this  is  seen  in  Fig«  29 


Figure  29 


In  terms  of  the  discrete  time  samples  we  have 

Xo(A-m)=  /S  ^  ^  ^vi  ( 


—  ccrvy^ 


T  g  sample  period 
a  s  time  constant 

The  function  was  truxjkated  at  n  §  20  for  values  of  k  ^  20« 
The  contribution  to  the  filtered  value  from  the  terms  multiplied 
by  e"^^  for  n  ^  20  is  negligible. 
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The  coherency  Is  a  statistical  parameterj  and  raist  be  repre= 
sented  as  such. 

Let? 


X  =  «1  -  «2 

A  2  -  ”1  -  “3 

^  3  s  H2  =  H3 


The  differences  of  the  unfiltered  values  were  grouped  in  a  frequency 
distribution  which  approximates  the  continuous  density  function  for 
large  sample  size.  Since  the  process  is  digital^  the  allowable  dif¬ 
ferences,  in  bits  ares  0,  ^  1,  ^  2,  ^  »«<■.  Therefore,  the  frequency 
distribution  is  not  a  smooth  cvirve.  Since  the  instrument  noise  is 
randomly  distributed,  the  envelope  of  the  distribution  is  a  normal 
density  function.  The  mean  was  removed  before  plotting,  and  the 
statistic  representing  the  standard  deviation  was  computed  according 


where s  N  =  the  ntuaber  of  samples  in  a  6.5  hour  record* 


Comparison  of  the  density  function  curve  with  an  ideal  signal 
input  to  that  of  magnetometer  input  gives  a  measure  of  the  noise 
produced  by  the  device. 

To  f’urther  define  the  coherence,  the  auto  and  cross  correla¬ 
tion  functions  were  computed  and  the  coherency  computed  as  defined s 
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where s 


and 


L 


L 


This  discrete  computational  form  is 


and 


where  s 

N  g  number  of  data  values 
n  3  number  of  lags 
i  -  n 


Since  the  functions  are  from  finite  time  series  can  only  be  made 

about  0.1  IJ.  Tho  6.5  hour  records  will  vary  in  form  deperdirsg  upon 
the  time  of  day  the  records  were  made.  Therefore,  the  correlation 


functions  will  not  be  the  same  shape  for  each  reeord|  however  the 
coherence  function  as  defined  will  be  the  samej  and  for  perfectly 
coherent  records  is  unity*  The  programs  are  included  in  Appendix 
III. 

It  should  be  pointed  out  that  when  perturbations  are  placed  in 
the  field,  the  basic  assumptions  regarding  correlation];  ioe*  a  sta^ 
tionary  random  process,  are  no  longer  applicable* 

Experimental  data  was  taken  to  determine  typical  signatures  for 
elementary  perturbations.  Fig®  12  shows  the  geometry  of  the  detector 
centered  coordinate  system® 


PERTURBATION  EXPERIMENT  GEOMETH? 
1- 


Figure  12 
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The  dipole  moment  was  orientated  parallel  to  the  earth's  field  and 
its  motion  was  in  the  x  direction  such  that  dx/dt  =  constant#  The 
approximation  r  1  is  valid# 

At  the  detectors 

H  g  ]j  Ho^  =  Hp^ 


Ands 


Hp 


where s 

M  s  dipole  moment 

r  s  function  of  time  and  the  geometry 
K  s  constant 

No  attempt  was  made  at  this  time  to  compare  the  measured  values  with 
the  theoretical  values ^  because  the  terrain  does  not  lend  itself 
keeping  L,  d,  and  v  constant.  The  results  are  representative  of 
typical  signatures  that  X'Jould  be  found  in  practice. 
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U.  Results 


Ij.l  Instrxanenta  tloiia 

An  ideal  signal  was  substituted  for  the  procession  signal  and 
set  of  records  were  made,  as  shown  in  Figs.  13  through  The 
graph  scales  were  left  unchanged  because  it  was  desired  to  show  the 
instrument  errors  on  the  same  scale  as  the  magnetic  field  scales. 

The  magnetic  field  vs.  time  graphs  were  not  included  because  the  in¬ 
put  signal  was  an  audio  oscillator  and  only  a  small  drift  puts  it 
off  the  graph  scale.  The  clock  was  300  kc/s  which  makes  1  bit  a 
0.36  /  . 

Figure  28  shows  the  precesion  signal  from  the  various  sensors. 
Figs.  26  and  27  are  typical  plots  made  from  the  analog  section  of  the 
instrument.  For  comparison,  the  Rb®7  magnetometer  output  was  adjust¬ 
ed  to  the  same  scale  and  plotted  on  the  same  graph.  The  Rb®^  magneto¬ 
meter  output  is  an  analog  of  frequency,  while  the  proton  magnetometer 
the  analog  is  of  period^  hence  the  reason  for  the  opposite  drift  on 
the  chart.  Figure  27  is  an  analog  record  of  the  differences.  The 
filter  was  a  single  section  RC  filter  with  time  constant  §  10  sec. 

On  the  miiLtiple  curve  plots  showing  coherence,  the  curves  are  identi¬ 
cal  and  cannot  be  distinguished  from  each  other.  These  curves  can 
be  compared  with  the  curves  obtained  from  the  sensor  with  no  local 
pertTirbations  in  the  field,  for  the  contribution  of  the  measuring 
instrument  noise  to  the  total  noise. 

Difficulties  with  both  instrumentation  and  computer  tape  have 
been  encountered.  The  use  of  new  computer  tape  for  data  storage  and 
scratch  tapes,  and  frequent  head  cleaning  of  the  instrumentation 
recorder  reduced  the  difficulties  to  a  minimumc 


29 


yigvtre  13 

X-SCfiLE  =  1.08Z+B3  miTS/INCH. 

Y-SCfiLE  -  Q.aSE+Oa  IMT3/'INCH. 

fiNDERSON  BOX  263 

first  difference  in  MRG  FIELI 


T  IN  HRS  H  GRMMR 


30 


Gi 


Differences  in  magnetic  intensity  vs. 
time  for  an  ideal  signal  input. 

DEL  1  s  Hi  -  H2 
DEL  2  =  Hi  -  H3 
DEL  3  r  H2  -  H3 

Filter  a  Z  transform  low  pass  with 
time  constant  of  20  minutes. 


<0 

CO 


? 


Figure  lit 

X-5Ca£  -  1.8ac+03  UHITB/'IHCH. 

Y-SCftL£  - 

^MDERSON  FILTER 

TLTEREDDIFFERENCE  in  MflG  FIELD  T  IN  HRS  H  GRMMfl 

31 


Figure 

}<-SC^LE  -  5^SE+e0  miTS/INCH. 

Y-SCfiLE  -  2.e8£-01  IMTS-'INCH. 

fiNDERSON  BOX  263 

density  function  X  IN  GflMMR  Y  IN  FREQ 


Normalized  density  function  for  the 
difference  Hi  -  H3  -  mean  (Hi  -  H3)>  ^ 

with  an  ideal  signal  input.  c& 

Kean  a  0.056  gamma 
Standard  deviation  r  0.72 
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Normalized  density  function  for  the 
difference  Hp  -  H3  -  mean  (H2  -  H3),  ^ 

with  an  ideal  signal  input.  ® 

Kean  3  —0.0055  gamma 
Standard  deviation  =  0,72 
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Cross  correlation  for  an  ideal 
signal  input. 

Record  length  r  6,5  hours 
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Figure  19 
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Coherence  for  an  ideal  signal  input 
Record  length  s  6, S  hours 
H]_2  =  Coherence  of  Hx  and  H2 
H13  =  Coherence  of  Hi  and  H3 
H23  =  Coherence  of  H2  and  H3 


m: 


^  esi  082  633  006  608  63^ 

Figure  20 

K-SCaE  -  :.08£-<-01  IHITS/^INCK 
'f-SCfLE  -  2.0e£-01  IMTS/JMCH. 

ANDERSON  BOX  263 

COHERENCE  FUNCTION  Y  IN  COH  X  IN  LAGS 

37 


-025 _ -620 _ -OJS  -010 _ -035 _ hgO _ 0^ _ 010 _ 015 _ ^ 


Differences  in  magnetic  intensity 
for  an  ideal  signal  input 
Record  length  =  1  hour 
Filter  »  2  point  average 
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Figure  21 


838 


38 


-ojfl _ -ess _ -e0i _ -092 _ ji:ii _  _ ggs _ ^ 


Difference  in  magnetic  intensity 
for  an  ideal  signal  input 
Record  length  s  1  hour 
Filter  =  Z  transform  low  pass 
with  time  constant  of  100  seconds. 
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Figure  22 
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Auto  correlation  for  an  ideal 
signal  input. 

Record  length  -  1  hour 
H]_  r  Correlation  (Hi  -  mean  Hi) 
H2  s  Correlation  (H2  -  mean  H2) 
H3  =  Correlation  (H2  -  mean  H3) 
Lag  -  A  t 
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Figure  23 
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Cross  correlation  for  an  ideal  signal  input. 
Record  length  3  1  hour 

Hi2  s  Correlation  (Hi  -  mean  Hi  and  H2  -  mean  H2) 

H13  m  Correlation  (Hi  -  mean  Hi  and  H3  -  mean  H3) 

H23  =  Correlation  (H2  -  mean  H2  and  H3  -  mean  H3) 
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Figure  2h 
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Coherence  for  an  ideal  signal  input 

Record  length  3  1  hour 

Hi2  -  Coherence  of  Hi  and  H2 

H13  a  Coherence  of  Hi  and  H3 

H23  =  Coherence  of  H2  and  H3 
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TYPICAL  ANALOG  PLOT  OF 
MAGNETIC  INTENSITY  VS.  TIME 
FILTER  .  RC  WITH  TIME  CONSTANT  OF  6  SECONDS 


FIGURE  26 


TYPICAL  ANALOG  PLOT  OF 
THE  DIFFERENCES  VS.  TIME 
FILTER  »  RC  WITH  TIME  CONSTANT  OF  6  SECONDS 
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FIGURE  27 
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li«2  Magnetic  field  moasvtreiaents 

Only  representative  results  of  the  magnetic  field  measuranients 
are  presented^  since  the  measurements  are  not  complete  at  this  time* 
Samples  of  plots  riade  with  the  clock  at  100  kc/s  are  included  to 
illustrate  the  effect  of  the  clock  frequency  on  the  results*  For 
the  associated  density  function  of  the  differences,  the  allowable 
values  are  0  ^  1*1,  ^  2*2,  •  ®  <>  in  gamma,  v;ith  the  clock  at  100 

kc/s.  The  expanded  difference  vs.  time  shows  clearly  the  ^  last 
bit  limitation*  Comparing  the  density  functions  with  the  plots  for 
the  clock  at  300  kc,  shov/s  that  most  of  the  noise  associated  wi.th 
the  measurement  is  due  to  the  last  bit  and  can  be  reduced  ty  increas¬ 
ing  the  clock  frequency*  The  perturbation  experiments  were  conducted 
with  the  clock  at  300  kc/s*  The  resuirs  of  the  perturoation  experi¬ 
ments  show  clearly  the  effect  of  the  digital  filtering,  and  the  im- 
poirbance  of  rime  constant  consideration  for  target  detection.  The 
oscillator  is  to  be  run  1  mc/s  for  further  experiments.  Also,  inapy 
of  the  errors  that  appear  in  the  data  were  the  direct  result  of  the 
instrumentation  tape.  Mew  tape  has  been  purchased  and  should  elim¬ 
inate  most  of  these  errors,  in  future  records. 

As  pointed  out  previously,  the  basic  assumptions  regarding  cor¬ 
relation  are  not  valid,  vrhen  the  noise  is  not  stationai-yj  however, 
since  correlation  detectors  can  be  instrumented  and  may  prove  useful 
in  the  detection  problem.  Some  of  the  eon'elation  function  plots  are 
included  to  illustrate  the  effect  of  local  pertui'bations  on  the  func¬ 
tions. 
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Conclusions 


For  target  detections  where  the  signals  under  consideration  are 
in  the  ^  gammas  or  less  ranges  it  will  be  difficult  to  detect  the 
target  •without  resorting  to  sophisticated  signal  processing*  In 
this  range  of  signalss  consideration  should  be  given  to  other  instru= 
mentsj  such  as  the  Rubidium  optically  pumped  luagne tome ter s  whose 
theoretical  sensitivi-ty  is  OoOOl  gammas  especially  where  the  detector 
is  stationary’s  orientation  is  not  a  problen*  It  can  be  demonstrated 
that  when  using  the  same  measurement  technique  as  for  the  proton 
magnetometers  the  instrument  and  processing  noise  will  be  an  order  of 
magnitude  less  than  the  theoretical  limits  imposed  by  the  sensiti- 
vity  of  the  instrument. 

For  signals  in  the  J  to  1^  gamuia  ranges  detection  could  be  ac" 
complished  with  analog  type  correlation  detectors.  These  could  be 
implemen'bed  with  operational  amplifiers. 

Signals  above  this  ran^e  could  be  detected  with  the  simple  dif- 
ferencing  techniques  and  the  associated  RcC,  filters. 

Increasing  the  clock  frequency  and  averaging  the  period  over 
more  than  1G2U  cycles  would  reduce  the  laeasm’ement  noise  to  a  smaller 
part  of  the  to’tal  noise.  The  target  signals  could  then  be  reduced 
by  a  factor  of  three^  with  no  clsange  in  the  Instrujr.en'bation. 

Full  utilization  of  the  inherently  narrow  line  width  of  proton 
magnetometer  cannot  be  achieved  using  the  standard  measuring  teeh= 
nique.  Other  techniques  have  been  suggested  by  Professor  Carl  I. 
Kennekenj  such  as  utilising  a  rapid  coherent  polarization  in  a  phase* 
lock  loop  arrangement*  Further  devolopiuent  work  needs  to  be  don^ 
to  explore  these  possibilitieso 
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•hour.  A  5  gamma  step  applied  to  at 
t  r  1  hour  and  removed  at  t  -  1  3/u  hours. 


m 


mi 


e82 


833 


50’i 


oas 


03Q 


oa^ 


-  1.03E+01  IMTS^IMCH. 
Y-5CflLE  »  Q.03E-01  UNITS/IHi^ 


ANDERSON  BOX  263 

COHERENCE  FUNCTION  Y  IN  COH  X  IN  LAGS 


1-36 


APPEI'IDIX  II 


16Q  FTOGRAMS  AllD 
OPERATOR  IHSTRUCTIOI'iS 


II 


APPENDIX  II 


160  PROGRAMS  AND  OPERATOR  INSTRUCTIONS 


Starting;  Addresses 


Function 


0000 

0002 

0004 

0006 

0010 

0100 


Rewind 

Read 

Search 

Write 

Write  identification  Block 
Main  Input  Program 


Error  Stops 


Cause 


ERROl 

ERR02 

ERR04 

ERR03 


773  of  current  input  in  error 
Memory  storage  not  complete 
End  of  storage  tape 
163  not  ready 


Addresses  of 
Important  Counters 


Information 


0057 

0045 

0046 


Number  of  blocks  v/ritten  on 
storage  tape 

Current  storage  location  in 
memory 

Number  of  zero  words  in  mem¬ 
ory 


Operator  Instructions 

1.  Load  program  and  make  the  necessary  cable  connections o 

2.  Set  Identifying  records  in  cells  3O3  to  370. 

3.  Write  identifying  block  (run  from  P  s  010) «  Stop  at 
P  s  100. 

4.  Set  cell  100  s  2200 

57  s  0 

5.  Start  Ampex  C«P,  100  and  run  from  P  g  100,  Stop  when 
tape  is  finished  or  ERR  stop. 
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0100 

0101 

0102 

0105 

0104 

0105 

0106 

0107 

0110 

0111 

0112 

0115 

0114 

0115 

0116 

0117 

0120 

0121 

0122 

0125 

0124 

0125 

0126 

0127 

0150 

0151 

0152 

0155 

0154 

0155 

0156 

0157 

0140 

0141 

0142 

0145 

0144 

0145 

0146 

0147 

0150 

0151 

0152 

0153 

0154 

0155 

0156 
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2200 

LDFOO 

1000 

1000 

4045 

STD45 

(ADDR) 

0400 

LDNOO 

4046 

STD46 

(GNTRl) 

0400 

LDNOO 

4150 

ST  1 50 

4151 

ST  1 51 

4152 

ST  1 52 

4155 

ST  1 53 

4154 

ST  1 54 

4155 

ST  1 55 

7500 

SXFOO 

0500 

0500 

7205 

INPO3 

0620 

0620 

6102 

NZF02 

0600 

0600 

1600 

LSFOO 

0607 

ADNO7 

6115 

NZFI5 

2150 

LDI50 

6015 

ZJFI5 

2151 

LDI5I 

6011 

ZJFll 

2152 

LDI52 

6007 

ZJF07 

2153 

LDI55 

6005 

ZJFO5 

2154 

LDI54 

6005 

ZJFO3 

2155 

LDI55 

6105 

NZFO5 

5446 

A0D46 

077Y 

SBN77 

6536 

NZB36 

0001 

EElROl 

2050 

LBD50 

4047 

STD47 

2147 

LDI47 

4145 

ST  145 

5445 

A0D45 

5447 

A0D47 

3456 

SBD56 

6505 

NZB05 

2200 

LDFOO 

4777 

4777 

II=.5 


0157 

0160 

0161 

0162 

0163 

0164 

0165 

0166 

0167 

0170 

0171 

0044 

0045 

0046 

0047 

0050 

0051 

0052 

0053 

0054 

0055 

0057 
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3445 

6653 

2044 

3445 

6003 

6657 

0002 

7101 

0172 

1000 

7700 


SBD45 

PJB53 

LDD44 

SBD45 

ZJF03 

PJB57 

EEIR02 

JFIOl 

WRITE 


HLTOO 


ADDITIONAL  CELLS  USED 


(PUT  IN  LIMIT) 
7600 


LIMIT 

ADDR 


CNTR  1 
CNTR  2 

BUFFER  STORAGE 
BUFFER  STORAGE 
BUFFER  STORAGE 
BUFFER  STORAGE 
BUFFER  STORAGE 
BUFFER  STORAGE 

NUI4BER  OF  MEMORY  DUMPS 
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0172 

0173 

0174 

0175 

0176 

0177 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

0210 

0211 

0212 

0213 

0214 

0215 

0216 

0217 

0220 

0221 

0222 

0223 

0163 

0164 

0165 

0166 

0167 

0170 

0171 
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5457 

A0D57 

2600 

LCFOO 

7776 

4060 

STD60 

7500 

EXFOO 

1141 

REQ  STAT 

7600 

INAOO 

0202 

LPN02 

6004 

ZJF04 

5460 

A0D60 

6506 

NZB06 

0003 

ERR03 

7500 

EXFOO 

2111 

2111 

7303 

0UT03 

7600 

TEEIM  OUT  (LI 

6102 

NZF02 

1000 

INIT  WRIT 

7500 

EXFOO 

1141 

REQ  STAT 

7600 

INAOO 

0240 

LPN40 

6103 

NZF03 

7101 

JFIOl 

0100 

00 

7700 

HLTOO 

+  1) 


DUI'^P  TO  PAPER  TAPE 


7500 

4104 

7303 

7600 

6102 

1000 

7700 
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DAILY  IDENTIFICATION 


USE  8  CELLS 

0030 

0070 

EXAMPLE 

3 

0031 

0064 

0032 

0032 

J 

0033 

0030 

A  month 

0034 

0060 

N 

0035 

0072 

6 

0036 

0062 

4  year 

0037 

0045 

CAR  RET 

11=7 
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S'p/JRT’ 


DAILY  RECORD  TO  M.T. 


i~p£> 
S  O 


S'PX 

looo 


0224 

0225 

0226 

0227 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0260 

0261 

0262 

0263 

0264 

0265 

0266 

0267 

0270 

0271 

0272 

0273 

0274 
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DAILY  RECORD  TO  MoTo 


2030 

LDD30 

4100 

ST  100 

1000 

1000 

2031 

LDD31 

4100 

ST  100 

1001 

1001 

2032 

LDD32 

4100 

ST  100 

1002 

1002 

2033 

LDD33 

4100 

STIOO 

1003 

1003 

2034 

LDD34 

4100 

STIOO 

1004 

1004 

2035 

LDD35 

4100 

STIOO 

1005 

1005 

2036 

LDD36 

4100 

STIOO 

1006 

1006 

2037 

LDD37 

4100 

STIOO 

1007 

1007 

2200 

LDFOO 

7700 

7700 

4100 

STIOO 

0100 

0100 

2200 

LDFOO 

1010 

1010 

4045 

STD45  (ADDR) 

0400 

LDNOO 

4145 

ST  1 45 

2044 

LDD44 

3445 

SBD45 

6003 

ZJF03 

5445 

A0D45  (ADDR) 

6506 

NZBO6 

7101 

JFIOl 

0172 

0172  (DU14P) 

7700 

NORl-IAL  STOP 

II-9 


REVJIND,  READ  AND 


rki  i.  jL^ 


Rmm).  RFAD  AND 
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Start 


Rewind  0275 

0400 

LDNOO 

0276 

4043 

STD43 

0277 

7500 

EXFOO 

0300 

1141 

REQ  STAT 

0301 

7600 

INAOO 

0302 

0202 

LPN02 

0303 

6504 

NZB04 

0304 

7500 

EXFOO 

0305 

1161 

REWIND 

0306 

7500 

EXFOO 

0307 

1141 

REQ  STAT 

0310 

7600 

INAOO 

0311 

0202 

LPN02 

0312 

6504 

NZB04 

Read  0313 

7700 

HLTOO 

Start  0314 

7500 

EXFOO 

0315 

1141 

REQ  STAT 

0316 

7600 

INAOC 

0317 

0202 

LPN02 

0320 

6504 

NZB04 

0321 

7500 

EXFOO 

0322 

2111 

EXFll  READ 

0323 

7203 

0UT03 

0324 

7600 

TERM  READ 

0325 

6102 

NZF02 

0326 

1000 

INIT  READ 

0327 

7500 

EXFOO 

0330 

1141 

REQ  STAT 

0331 

7600 

INAOO 

0332 

0240 

LPN40 

0333 

6002 

ZJF02 

0334 

0004 

ERR04  (END 

0335 

2043 

LDD43  TA; 

0336 

0701 

SBNOl 

0337 

6004 

ZJF04 

0340 

6525 

NZB25 

Search  0341 

0401 

LDNOl 

Start  0342 

4043 

STD43  (  SEA] 

0343 

2100 

LDIOO 

0344 

1000 

1000 

0345 

3430 

SBD30 

0346 

6532 

NZB32 

0347 

2100 

LDIOO 

0350 

1001 

1001 

0351 

3431 

SBD31 

0352 

6536 

NZB36 

0353 

2100 

LDIOO 

11-12 
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0354 

1002 

1002 

0355 

3432 

SBD32 

0356 

6542 

NZB42 

0357 

2100 

LDIOO 

0360 

1003 

1003 

0361 

3433 

SBD33 

0362 

6546 

NZB46 

0363 

2100 

LDIOO 

0364 

1004 

1004 

0365 

3434 

SBD34 

0366 

6552 

N2B52 

0367 

2100 

LDIOO 

0370 

1005 

1005 

0371 

3435 

SBD35 

0372 

6556 

NZB56 

0373 

2100 

LDIOO 

0374 

1006 

1006 

0375 

3436 

SBD36 

0376 

6114 

NZF14 

0377 

2100 

LDIOO 

0400 

1007 

1007 

0401 

3437 

SBD37 

0402 

6110 

NZFIO 

0403 

2100 

LDIOO 

0404 

1010 

1010 

0405 

6105 

NZFO5 

0406 

0400 

LDNOO 

0407 

4043 

STD43 

0410 

7101 

JFIOl 

0411 

0077 

0077 

0412 

7101 

JFIOl 

0413 

0314 

TEST  TAPE  22  BLOCKS 

DATA 

M  *  4714510 

LIMIT 

CELL 

44  = 

7600 

ADDR 

CELL 

45 

-  0013  40518 

CNTR  2 

CELL 

21 

N2  Ml 

Ml 

CELL 

22 

N2 

CELL 

23 

CNTR  1 

CELL 

20 

0526 

7700 

HLTOO 

0527 

0400 

LDNOO 

0530 

4021 

STD21 

0531 

2200 

LDFOO 

0532 

1000 

1000 

11-13 


0533 

0534 

0535 

0536 

0537 

0540 

0541 

0542 

0543 

0544 

0545 

0546 

0547 

0550 

0551 

0552 

0553 

0554 

0555 

0556 

0557 

0560 

0561 

0562 

0563 

0564 

0565 

0566 

0567 

0570 

0571 

0572 

0573 
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4045 

0413 

4023 

0400 

4020 

2200 

4051 

4022 

2021 

0726 

6417 

2023 

4145 

5445 

2022 

4145 

5445 

5422 

5420 
2200 
4777 
3445 

6207 

2044 

3445 

6104 

5421 
7101 
0172 
2020 
0777 
6434 

6525 


STD45 

LDN13 

STD23 

LDNOO 

STD20 

LDFOO 

4051 
STD22 
LDD21 
PPN26 
ZJBI7 
LDD23 
STI45 
A0D45 
LDD22 
ST  1 45 
A0D45 
A0D22 
A0D20 
LDFOO 

4777 

SBD45 

PJFO7 

LDD44 

SBD45 

NZF04 

A0D21 

JFIOl 

0172 

LDD20 

SBN77 

ZJB34 

NZB25 
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WRITE  TEST  TAPE 


22  BLOCKS  DATA 
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FORTRAN  PROGRAMS 

ANALYSIS  PROGRIIi 

SATELLITE  DISPUY  PRXRAM 


UIMtNSlCJN  NDATA  (664)  ,  DATE(  00  ),  H(  17  28) ,  HONE  (  900  ,  IT  I  TLE  (  1  2)  , 
1HTW0(900) ,HTHRE (900) , T ( 9CC 1 , 8  IN  1 ( 10 0 ) » E I N2 (  1 00 )  , B I N3 (  1 00 ) , X ( 1 00 ) 
no  5  J  =  1,14 
READ  5,DATi-(J) 

5  E0RNAT(C16) 

COr>!(NA.SK  =  000000007777777  j'B) 

ENA(KDATA) , INA( 1 ) . 

STA( INI  D  ,  INA(664) . 

STA (  ITlRN) ,LNI ( 0) . 

6  CALL  READ(:.>,  INI  T,  ITERK,  1  ) 

DO  70  1=3,.- 

SLSl (7) 

7  AOATE  =  DATE ( I ) 

BDATE  =  DATE ( 1+ 1 ) 

LDA(  ADATE)  ,SL!R(  NDATA  +1  )  ,  AJPl  (6)  , 

LDA(BDAT6) ,SU8( NDATA  +2), AJPl (6). 

NBLOC  =  0 

30  CALL  READ(o, INI T, ITERM,  1 ) 

NZRC  =  0 
DO  81  J=1 ,864 
IF(NDATA(J))  82,82,81 
82  NZRO  =  NZRO+1 
81  CONTINUE 

NO  =  NDATA! 5) 

LDA(ND) ,AJP0( 10 ). 

K  =  C 

DO  40  J=1 ,864 
K  =  K+1 

INDATA  =  NDATA! J) 

LDA(  INDATA)  ,  LDG  ( MASK  )  ,  ARS  {  4' 4  )  ,STL(INLIST). 

H(K)=INLIST 
K  =  K+1 

LDA( INDATA) , LDQ (MASK ),STL(1NLIST). 

SLSl  (40  . 

40  H{K)=INL1ST 
AVt=C. 

DO  50  J  =1 , 1728 
H(J)  =  H(J)/3. 

H(  J )  =  1024./H( J ) *2348400 .0 
50  AVE  =  H(J)  +  AVE 


PRINT  26,  AVE 
C  ZAPP  FILTER 

NERR  =  0 

AVE  =  0. 

HBAR  =  50950.0 
DO  51  J=1 ,1728 
DEV  =  A8SF(H( J) -HBAR) 

IFIDEV-200.)  51,52,52 
JCGUNT  =  J 

52  NERR  =  NERR  +  1 
IFIJ-3)  53,53,54 

53  H( J)  =  H( J+3) 

GO  TO  51 

54  H(J)  =  H(J-3) 

51  AVE  =  AVE  +  H(J) 

AVE  =  AVE/1720. 

PRINT  26,  AVE 
PRINT  97, NERR, NZRO, JCOUNT 
NBLOC  =  NBLOC  +  1 
WRITE  TAPE  2,  H 
GO  TO  30 
10  REWIND  2 

C  THIS  SECTION  COMPUTES  H  VS  TIME 
C  18  PT  AVE  N  IS  TOTAL  NUMBER  CF  POINTS  PER  DAY 
TIME  =  -  6.5/704. 

DO  60  L=l, NBLOC 
READ  TAPE  2,H 
DO  62  J=1 ,32 
HIBAR  =  0. 

H2BAR  =  0. 

H3RAR  =  0. 

DO  61  K=1 ,16 

M  =  3*K  -  2  +  (  J  -  1  )*5  4  III-l 

HIBAR  =  H(M)  +  HIBAR 
H2bAR  =  H(M+1)  +  H2BAR 

61  H3BAR  =  H(M+2)  +  H3BAR 

N  =  J  +  (L-1 )*32 

HONE(N)  =  Hi  BAR/ia. 

HTWC(N)  =  H2  BAR/ 1-8. 

HTHRE(N)  =  H3BAR/18. 

TIME  =  TIME  +  6.5/704. 

62  T(N)  =  TIME 

_ _  r _ _ _ _  _ _ _ _ _ 


RE^^IND  2 
PR  1. N'T  2  7,  N 
PRINT  27, M 
NPTS=N 

C  MEAN  VALUE  OF  DIFFERENCES 
JERR  =  0 

I  ERR  =  0  APPENDIX  III 

NERR  =  0  - 

DELI  =  0. 
nEL2  =  0. 

DEL3  =  0. 

DO  64  J  =  1,NBLCC 
READ  TAPE  2,h 
DO  6U  L  =  1,576 
M  =  3*L-2 

DEL  =  H(M  +  1)  -  HIM) 

IF! ABSF (DEL)-  UO. )  65,66,66 

65  DELI  =  DELI  +  DEL 
GO  TO  67 

66  NERR  =  NERR  +  1 

67  DEL  =  HIM  +  1 )  -HIM  +  2  ) 

IFIABSFIDEL)  -  40.)  68,69,69 

68  0EL2  =  DEL2  +  DEL 
GO  TO  71 

69  I  ERR  =  lERR  4  1 

71  DEL  =  HIM)  -  HIM  +  2) 

IFIABSFIDEL)  -  40.)  72,73,73 

72  DEL3  =  DEL:i  4  DEL 
GO  TO  64 

73  JERR  =  JERR  4  1 
64  CONTINUE 

REWIND  2 

PRINT  97, NLRR, I  ERR, JERR 
97  FORMAT! 3120 
ABLGC  =  NBLOC 
DELI  =  DELI / (576. *AELOC  ) 

DEL2  =  DEL2/ 1576. «A8L0C ) 

DEL3  =  DELS/ 1576. »ABLOC ) 

PRINT  49, DELI ,DEL2, DELS 
C  DENSITY  FUNCTION 
STDEV1=0. 

ST0EV2=0. 

STDEV3=0. 

NERR  =  0 
lERR  =  0 
JERR  =  0 
XL=-25. 

DO  45  K=l,100 
B  I N 1  I K ) =C . 

BIN2(K)=0. 

BIN3I K) =0. 

X(K)=XL 
45  XL=  XL40.5 

DO  74  J=l, NBLOC 
READ  TAPE  2,h 
DO  74  L  =  1,576 
M  =  3»L-2 

STDEV1=  IH(M4l )-HIM)-CEL1 )*e2  4STDEV1 
STDEV2=  (H(M4l)-HIM42)-CEL  )*42  4  STnEV2 
STU£V3=  (H(M)-HIM42)-DEL5)»*2  4STDEV3 
DEL  =  H(M4l )  -  HIM) 

IF! ABSF I  DEL-DEL l)-25. ) 7 5, 75, 76 

76  NERR  =  NERR  4  1 

75  DEL  =  HIM41 )  -HIM42) 

IF! ABSF I  DEL-DEL 2) -25. ) 7 7, 77, 78 
78  IFRR  =  lERR  4  1 

77  DEL  =  HIM)  -  H(M42) 

IF! ABSF I  DEL-DEL  3) -25. ) 2 9,  29,80 
80  JERR  =  JERR  4  1 
29  DO  31  K=1 ,3 

IF(K-2)32,33,34 

32  DEL=H(M4l )-H(M) -DELI 
GO  TO  35 

33  DEL=HIM4l )-H(M42)-DEL2 
GO  TO  35 

34  DEL=H(N)-HIM42)-DEL3 

35  XL=  -25. 

DO  36  N=1 , 1 00 

IFIDEL-XL)  37,37,38  III-2 

38  IF(DEL-(XL40.5) )39,39,37 

39  IF(K-2)  42,43,44 

42  PINl (N)=BIN1 (N)4l . 

GO  TO  31 

43  D:N2(N)  =  BIN2(N)4l. 

GO  TO  31 
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O  O  U  U  K  I  1  l\  u  t 
31  CCNTINbE 
7U  COMINUh 

STbEVl  =SQRTF(STDEV1/{576.=ABL0C) ) 

STDEV2  =SQRTF(ST0EV2/(576.«ARL0(:)  ) 

ST1)EV3  =SQRTF(STDEV5/ (  57d.«  A6L0C)  ) 

PRINT  Ug.STDEV  1  ,STDEV2,  STD£;V3 
PRINT  97,NtRR, lERR, JERR 

remind  2 

NORMALIZED  DENSITY  FUNCTION 
AMAX=BIN1 ( 1 ) 

BMAX=BIN2{ 1 ) 

CMAX  =  P,IN3{  1  ) 

no  200  J=1 , 1 00 

IF ( AM AX-BIN  1  ( J)  )201 ,202,20 

201  AMAX=RIN1{J) 

202  IF! BMAX-eiN2( J)  )203,204,20L 

(-) +  +  t  +  *  +  J  +  +  +  t  +  S  +  t ♦  ++=!=  +  //++  +  +++ +  +  +  +  5 

20U  IF (CMAX-BIN3 { J) 1205,200,200 
205  CMAX=BIN3(J) 

200  CONTI  NUN 


// 


PRINT  U9, AMAX, BMAX,CMAX 
DO  206  J=1, 100 
BIN1 ( J) =BIN1 ( J) /AMAX 
BIN2( J)=DIN2( JI/BMAX 
206  BIN3( J)=BIN3( J) /CMAX 
DO  46  K  =  1 , 12 
46  ITITLE(K)=8H 

ITITLEI 1 )=CHANO£RSCN 
ITITLE(2)=8H  BOX  2o3 
ITITLE{7)=£H  DENSITY 
ITITL£(8)=aH  FUNCTIO 
IT  ITLE(9)=6HN  X  in  G 
ITITLEI 10)=8HAMMA  Y 
ITI TLEII 1 )=8HIN  FREC 
LA6A=4H  pi 
LAR8=4H  P2 
LAbC=4H  P3 

CALL  DRAW(100,X,BIN1 ,C,C,LAbA, ITI TLE, 0,0, 0,4, 0,2, 8, 8,0, LA  ST) 
CALL  0RAW(|00,X,8IN2,0,C,LAfiB, I T ITL E , 0 , 0 , 0 , 4 , 0 , 2 , 8 , 8 , 0 , LA S T ) 
CALL  DRAW! 100,X,BIN3,0, C,LA8C, I  TITLE, 0,0, 0,4, 0,2, 8, 8,0, LA  ST) 
506  N  =  N PTS 

H0NEAV=0. 

HTWCAV=0. 


HTHREAV=0. 

DO  48  L=1,N 
HONEAV=HONeAV+HCNE{ L ) 
HTWCAV=HTWf;AV+H.TWC(L  ) 

48  HTHREAV=HTHREAV+HTHRE(L ) 

SV04  =  N 

HONEAV  =  H0NEAV/SV04 

HTWCAV  =  HTWOAV/SVOU 

HTHREAV  =  HTHREAV/SV04 

PRINT  49, HONEAV, HTWOAV, HTHREAV 

49  FORMAT! 3E20. 7) 

PRINT  27,NDLCC 

27  FORMAT! 17) 

HONEAV=  IHONEAV  +HTWOAV  +HTHREAV)/3. 
PRINT  26, HONEAV 
26  FORMAT! E20. 7) 

WRITE  TAPE  2,  HONE 
WRITE  TAPE  2,  HTWO 
WRITE  TAPE  2,  HTHRE 
REWIND  2 
DO  90  L  =  1,N 

HONEIL)  =  HONEIL)  -  HONEAV 
HTWCIL)  =  HTWOIL)  -  HONEAV 
90  HTHREIL)  =  HTHREIL)  -  HCNEAV 
DO  1  L=1,12 
1  ITITLEIL)  =  8H 

ITITLEI  1  )=8HANDERS0N 
ITITLE!2)=8H  BOX  263 
ITITLE!7)=3H  EARTHS 
ITITLE!8)=8HMAGNETIC 
ITITLE!9)=SH  FIELD  V 
IT  ITLE! 10)  =  8HS  TIME 
ITITLEI 1 1 )=8HT  IN  HRS 

ITITLEI  12)=8H  H  GAMMA  III-3 

LABA=4H  HI 
LABB=4H  H2 
LABC=4H  H3 


C  Z 


transform 

1.  I.ILIA/inn 


LAST  ) 
LAST  ) 
LAST  ) 


IF{J-20)  302,303,303 
502  M=J 

GO  TO  soil 

303  M=;o 

304  00  301  K=1,M  APPENDIX  III 

L=  K-1  - 


DL=L 

ASUI^  =  tXPFI -AI  *PL  )+ASUP 
BSUM  =1:  XPF(- AT  *61)  IP  HONE  (J-D  +  BSbK 
CSUK  =EXPF(-AT*eL) «HTWO(J-L )  +  CSUM 
301  OSUP  =  EXPr  ( -A1  «DL)  *HTHRE(  J-D+OSIM 
HONE(  J)=BSUM/ASUFi 
HTWO{J)=  CSUN'/ASUM 
300  HTHRE(J)=  DSUM/ASUM 
ITITLE(2)=BH  FILTER 

CALL  DRAW(N,  T,HCN!E,  1,0,LABA,ITITL£,  l.C,10.0,5,1,2,2,8,1C,0,LAST) 
CALL  ORAWIN, T,HTW0,2, 0, LABB,ITITLE,  1.0,  10. 0,5, 1,2, 2, 8, 10,0, LA ST) 
CALL  DRAW (N, T,HTHRE, 3,0 ,LABC, ITITLE, 1.0, 10., 5, 1,2, 2, 8, 10,0, LAST) 
C  DIFFERENCES  OF  FILTERED  VALVES 
DO  306  J  =  1 ,N 
H( 4)=  HONE( J ) 

HONE( J)=  H( J)-  HTWO( J) 

STORE  =  HTWO(J) 

HTWO(J)=  H{ J )-HTHRE{ J ) 

306  HTHRL(J)=  STCRE-  HTHREIJ) 

ITITLE(7)=6HFILTERED 
ITITLEt  8)=eH0IFFEREN 
TTITLE{9)=8HCE  IN  NA 
ITITLEI 10)=eHG  FIELD 
LAljA=4HD£Ll 
LABB=4HDEL2 


LABC=4HDEL3 

CALL  ORAW(N,T,HCNE, 1 , 0, L AB A , I T 11 LE,  1 . C , 5 . 00 , 5 , 1 , 2 , 2 , 3 , 1 0 , 0 , L A ST ) 
CALL  DRAW(N,  T  ,  HTUO,  2 , 0  ,  LABT! , ITITLE,  1.0,5.00,5,  1,2, 2, 8, 10,0,  LA  ST) 
CALL  DRAW(N,T,HTHRt,3,0,LA6C,  IT  I  TLE , 1 . 0 , 5 . 0 , 5 , 1 , 2 , 2 , 8 , 1 0 , 0 , LA ST ) 
C  THIS  SECTION  COMPUTES  SIMPLE  Di FF 
READ  TAPE  2,  HONE 
READ  TAPE  2,  HTKO  . 

READ  TAPE  2,  HTHRE 
REWIND  2 
DO  15  L=l,900 
H(L)=HONE{L) 

HONE(L) =H(L)-HTWC(L ) 

STORE=HTWO( L ) 

HTKO(L)=H{L)-HrHRE(L) 

15  HTHRE(L)=STORE-HTHRE(L) 

DO  16  L=1  , 12 

16  ,  ITITLE! L)=8h 

ITITLE!  1  )=eHANDERSQN 
ITITLE!2)=eH  BOX  263 
ITITLE!7)=8HFIRST  DI 
lTITLE!8)=eHFFERENCE 
ITITLE!9)=8H  IN  MAG 
ITITLE!  10)  =  8HFIELD 
ITITLE! 1 1 )=SHT  IN  HRS 
ITITLE! 12)=8H  H  GAMMA 
LABA=4HDEL1 


LABB=4HDEL2 

LABC=4HDEL3 

CALL  DRAW!N,T,HONE, 1 , 0, LAB A , ITITLE,  1.0,5. 00, 5, 1,2, 2, 8, 10,0, LAST) 
CALL  DRAW !N, T,HTWO, 2,0, L ABB , ITITLE,  1.0,5.00,5, 1,2, 2, 8, 10,0, LAST) 
CALL  DRAW ! N, T, HTHRE, 3, 0,L ABC, IT  I TLE , 1 . 0 , 5 . C , 5, 1 , 2 , 2 , 8 , 1  0 , 0 , L AST ) 
AUTOCORRELATION  FUNCTION  75LAGS 
NLAGS=75 
MPTS=  N-NLAGS 
PTS=  MPTS 


9 


1  1 


XL  =  0. 

DO  9  J=  1  ,  100 
B I N 1 ! J ) =0 . 
eiN2! J)=0. 

8IN3! J)=0. 

X!  J)=XL 
XL  =  XL+1. 

READ  TAPE  2, HONE 
READ  TAPE  2,HTW0 
READ  TAPE  2, HTHRE 
REWIND  2 
DO  11  J=1,N 

HONEIJ)  =  HONEIJ)-  HONEAV 
HTWOvJ)  =  HTWC!J)  -  HTWCAV 
HTHRE! J)=HTHRE! J)-HTHREAV 
DO  12  J=1,NLAGS 
DO  14  K=1,MPTS 
L=K-1 


r>  T  •>  >  *«  '  «  * 


r.tA  V  «-l  lO  •>  l.CL  ^  1  » 


O  T  t  ^  l_% 


tili\JtJ)=  blN^(J)  +  HTWO(K)<HTWO(L+J  ) 

14  f3Ii\3(  J)=RIN?(  J)  +  HTHRt{K)^HTHRL(L+J  ) 

BIN1 ( J)=BIK1 ( J) /PTS 
BI,\2(J)=  SIN2(J)/PTS 
12  BIN5(J)=  BIN3(J)/PTS 
C  SAVE  AUTOCGRRELATION  FUNCTION 
DO  105  J=l,100 

H(J)=eiN1(J)  APPENDIX  III 

H( J+100)=BIN2( J) 

105  H( J+200)=BIN3( J) 

00  100  J=l, 12 

100  ITITL£(J)=FH 
ITITLEI 1 )=ChANOERSCN 
ITITLE(2)=oH  BOX  263 
ITITLE(7»=:8HAUTCC0RR 
IT1TLE(8)=PHELAT10N 
ITITLEI 9) =3HFUNCTI0N 
ITITLEI 10)=0H  Y  IN  PR 
ITITLEI 1 1 )=8H00UCTS  X 
ITITLEI 12)=8H  IN  LAGS 
LABA=4H  HI 

LABD=4H  H2 
LABC=4H  Hi 

CALL  DRAW(NLAGS,X,BIN1t  1,0,LABA, I  TITLE, 0,0, 0,0, 0,0, 8, 8,0, LAST) 
CALL  DRAW(KLAGS,X,BIN2, 2,0,LABB, ITITLE,C,0,0,0,0,0,8,8,C,LAST) 
CALL  DRAWINLAGS ,X, BIN3, 3,0,LABC, I TI T L E , 0 , 0 , 0 , 0 , 0 , 0 , 8 , 8 , 0 , L AST  ) 
C  CROSS  CORRELATION  FUNCTION  75  LAGS 
DO  101  J=1,100 
BIN1 ( J)=0. 

BIN2( J)=0, 

101  PIN3(J)=0. 

DO  102  J=1,NLAGS 
DO  103  K=1,MPTS 
L=  K-1 

BIN1 ( J)=PIN1 ( J )+  HTWOIK )*H0N5(L+J ) 

BIN2 ( J) =8IN2 ( J ) +HTWO( K) «HTHRE(L+J) 

103  BIN3( J)=BIN3(J)+  HONE ( K  ) «HTHRE ( L  + J) 

BIN1 ( J)=BIN1 ( J)/PTS 

BIN2(J)=  BIN2(J)/PTS 

102  BIN5(J)=  B1N3(J)/PTS 
DO  104  J=l,12 

104  ITITLE{J)=SH 
ITITLEI 1 )=8HANDERS0N 
1TITLE(2)  =8H  BOX  263 
ITITLE(7)  =SHCRGSCaRR 
ITITLE(8)  =8HELATIGN 
ITITL£(9)  =OHFUNCTION 
I  TITLE!  10)=eH  Y  IN  PR 
ITITLE! n )=8H0CUCTS  X 
ITITLE(12)=eF  IN  LAGS 
LABA  =  4H  H21 

LABB  =  4H  H23 
LABC  =  4H  HI  3 

CALL  DRAWINLAGS ,X, 8  INI,  1,0, LABA, ITI TL E , 0 , 0 , 0 , 0 , 0 , 0 , 8 , 8 , 0 , L AST  ) 
CALL  DRAWINLAGS ,X,BIN2, 2,0, LABB, ITI TLE, 0 , 0 , 0 , 0 , 0 , 0 , 8 , 8 , 0 , L AST ) 
CALL  DRAWINLAGS,X,BIN3, 3,0, LABC, ITI TL E , 0 , 0 , 0 , 0 , 0 , 0 , 8 , 8 , 0 , L AST ) 
C  COHERENCY  COMPUTATION  75  LAGS 


107 

108 


109 

no 


1  1  1 
106 


DO  106  J  =  1,NLAGS 

BINl IJ)  =  IBIN1IJ)«*2)/IH{J  +  100)*HIJ)  ) 

IFIBINIIJ))  107,108,108 
BINI I J)=-BIN1 IJ) 

BINl I J) =  SORTFI BIN! IJ ) ) 

BIN2IJ)  =  I  BIN2IJ  )«*f2)/I  HIJ+i00)»HIJ+200)  ) 

IFI8IN2IJ))  109,110,110 
BIN2 ( J) =-BIN2 I J ) 

DIN2IJ)=  SORTFI BIN2IJ ) ) 

BIN3IJ)=I BIN3IJ ) ee2) / (HI J )^H(J+200) ) 

IFIBINSIJ))  111,106,106 
BIN3(J)=-EIN3IJ) 

BIN3(J)=  SQRTFI BIN3(J ) ) 

ITITLE(7)=eHC0H£RENC 
ITI TLEI 8) =8HE 
ITITLEI 10)=8H  Y  IN  CO 
ITITLEI 1 1 )=8hH  X  IN  L 
ITITLEI 12 )=8HAGS 

CALL  DR AW (NL AGS, X, BINl,  1,0, LABA, ITITLE,0,.2,0,0,0,0,8,8,0,LAST) 
CALL  DRAW(NLAGS,X,8IN2, 2,0, LABB, ITITLE,0,.2,0,0,0,0,8,8,0,LAST) 
CALL  DRAWINLAGS ,X,B IN3, 3, 0, LABC, I  title, 0, .2, 0,0, 0,0, 8, 8,0, LA  ST) 
CONTINUE 

REWIND  3  III-5 

STOP 


END 

MACHINE  READ! 1 ARG,2ARG, 3ARG,4ARG) 


TYPEWRITER  MESSAGE  CCOE  TABLE 


/irroivi-ixA  XXX 


C0N(T6  =4504041 22050220AP, 
«  TB  =121203 ! 2C0000C0Ce, 

*  110=1550 121 4C125042CP, 
CON (T12=450404 1220 30220 4E» 

*  T1U=141 1 2OOOOOOOOOOC0, 
140=4504040  1501520041', 

«  T42=C4140622141630016, 


T7  =11200613010504208, 
T9  =45040412203022048, 
T 11 =1 21 2031 2C00CC00CB ) 
T 13=20062204032604268, 
T15=4504C45112i40120B, 
T4 1 =5406 140 10406030 18, 
T 43 =20224200000000008) 


ft  TAPE  UNIT  ASSIGNMENT  TABLE 

ft  ASSIGNMENTS  FOR  CHANNEL  3/4  LIBRARY  TAPE 


CONIKO  =0 
ft  K3  =3320308 

ft  K6  =5520208 


K1  =3320^00, 
K4  =3320408, 
K7  =5520508, 


K2  =3320206, 
K5  =5520108, 
K8  =55^0408) 


ft  ASSORTED  CONSTANTS 

C0N{R4  =73737373737373758,  R 1 6=00700000000000008 ) 


ft  TAPE 

READING  ROUTINE 

IX 

SLJ(K) 

SLJ( L+5) 

.  EXIT/ENTRY 

1  ARG 

ZRU(O) 

ZRO(O) 

.  TAPE  UNIT  NUMBER  ARG. 

2ARG 

ZRO(O) 

ZRO(O) 

.  INITIAL  ADDRESS 

3ARG 

ZRO(O) 

ZRO( 0) 

.  TERMINAL  ADDRESS 

4ARG 

ZRO(O) 

ZRC( 0) 

.  MC05,  1  =  BINARY,  2  =  BC 

LDA7 ( 1 ARG) 

INA(-nB) 

SIUI ( 19X) 

AJP3( L+3 ) 

I  CHECK  FOR  TU  ASSIGNMENT 

ENA  (  T40  ) 

SLJ4(  G) 

.  GO  PRINT  ERROR 

SLJ(  P) 

ZRO( 0) 

INA{ IIB) 

STA(KO ) 

LILl  (KO)- 

LDA1 (KO ) 

ADD7(4ARG) 

SAU( 4X ) 

I  SET  TU  CONTROL  CODES 

SCL(777B) 

SAL( 4X ) , 

.  INITIALIZE 

SAL ( 5X) 

SAL( 1 IX) 

INA(3) 

SAU( BX  ) 

INA(2) 

SAU( 9X  ) 

INA(  1  ) 

SAu ( n  X ) 

INA{  1  ) 

SAU (7X  ) 

LRS (IS) 

ENA(  0) 

LLS(3) 

SAL( 2X ) 

SAL{3X) 

SAU ( 12X ) 

LLS(39) 

LDQ( R1 Z ) 

STA( KO) 

LDA( 5X ) 

SSU{ KO) 

STA( 5X ) 

LDA7(2ARG) 

SAU( 5X ) 

LDA(6X) 

SSU(K0 ) 

STA{6X) 

ENI ( 0) 

2X 

LDA7(3ARG) 

SAL( K) 

SAU( 6X) 

SIL2(  i9X) 

3X 

ENI2(4) 

L  I L  1  (  N  ) 

4X 

EXF (N) 

EXF7(N  ) 

I  SELECT  UNIT 

5X 

EXF ( N) 

EXF7(N ) 

.  ACTIVATE  CHANNEL 

6X 

EXF(N) 

LDA7(2ARG) 

SUh ( R4) 

LLS(40  > 

INI  1 (-1 ) 

LDAl ( Z  ) 

INI  1  ( 1  ) 

SUB(RA  ) 

7X 

EXF7(N) 

SLJ ( 13X) 

.*  SENSE  ENDFILE 

8X 

EXF7(N) 

SLJ( L  +  4  ) 

.  SENSE  PARITY  ERROR 

9X 

EXr7(N) 

SLJ(  lOX) 

.  SENSE  LENGTH  ERROR 

0JP1 ( 19X) 

AJP( 3X  ) 

.  SENSE  FOR  PAD  RECORD 

SLJ{ 19X) 

ZRC(O) 

.  EXIT 

J 

IJP2{ L+2) 

SLJ( 14X) 

.  SENSE  FIVE  PARITY  ERRORS 

lOX 

IJP2( L+  1  ) 

SLJ (  1  3X ) 

,  SENSE  FIVE  LENGTH  ERRORS 

AJP (3X) 

CJP( 3X  ) 

.  SENSE  FOR  BAD  RECORD 

1  1  X 

EXF (N) 

EXF7(N) 

.  BACKSPACE 

12X 

L  I  L  T  (  N ) 

SLJ ( 4X  ) 

.  RETURN  TO  REREAD 

13X 

£NA(T6 ) 

SLJ (  17X  > 

-  SET  TO  INDICATE  TROUBLE 

14X 

ENA{ T9) 

SLJ (  17X ) 

15X 

ENA ( T12  ) 

SL J( IBX) 

17X 

SLJ4(  G) 

ZRO(O)  III-6  .  TROUBLE  EXIT 

AJP ( 1 9X ) 

AJP3(5X) 

-  .  SENSE  ACTION  TO  TAKE 

SLJ  ( nx ) 

ZRO(O) 

18X 

SLJ4.;  g; 

ZRn( 0) 

I  TROUBLE  EXIT 

AJP ( 3X) 

AJP3(3X) 

.  SENSE  ACTION  TO  TAKE 

SLJ(llX) 

ZRO(O) 

1  9X 

ENI 1 (N) 

EN  I  2 ( N  ) 

.*  retStore  indexes 

SLJ (IX) 

ZRO( 0) 

PROGRAM  PROMAG  APPENDI;C  III 

DIMENSION  NDATA(864) ,DATE(  100) ,H( 1 7 2 8 ) , HON E ( 900 ) , I  TITLE (12), 
1HTW0(900) ,HTHRE(900) ,T( 900) 

DO  5  J  =  1,14 
READ  5,DATE(J) 

5  F0RMAT(016) 

CONIMASK  =  0000000077777777B) 

ENA(NDATA) , INA( 1 ) . 

STA(INIT).INA(864). 

STAI ITERM) ,£NI (0) . 

6  CALL  READ(.^,  INIT,  ITERM,  1 ) 

00  70  I  =  1,11,2 

ADATE  =  DATE ( I  ) 

BOATE  =  DATE (1  +  1  ) 

LDA( ADATE) ,SUB(NDATA  +1  ),AJP1 (6), 

LDA(BDATE) ,SUB(NDATA  +2),AJP1 (6). 

NBLOC  =  0 

30  CALL  READ(3, INIT, ITERM,  1 ) 

ND  =  NDATA(3) 

LDA( ND)  ,AJPO ( 10  )  . 

K=0 

DO  40  J=1,864 
K=K+1 

INDATA  =  NDATA(J) 

LDA( INDATA) , LDQ ( MASK ) , ARS ( 24 ) , STL ( I NLI ST ) . 

H(K)=INLIST 

K=K+1 

LDA( INDATA) , LDQ (MASK) , S TL ( I NL 1ST ) . 

SLSl  (40)  . 

40  H(K)=INLIST 
AVE=0. 

DO  50  J  =1,1 728 

H( J)=( 1024./H( J) ) *2348400.0 

50  AVE=H(J)+AVE 
AVE=AVE/1728. 

PRINT  26,  AVE 

C  ZAPP  FILTER 

NERR  =  0 

AVE  =  0. 

HBAR  =  50950.0 
DO  51  J=l, 1728 
DEV  =  ABSF(H( J)-HBAR) 

IF(0EV-200.)  51,52,52 

52  NERR  =  NERR  +  1 
IF(J-3)  53,53,54 

53  H(J)  =  H(J+3) 

GO  TO  51 

54  H(J)  =  H(J-3) 

51  AVE  =  AVE  +  H(J) 

AVE  =  AVE/1728. 

PRINT  26,  AVE 
PRINT  27,  NERR 
NBLOC  =  NBLOC  +1 
WRITE  TAPE  2,  H 
GO  TO  30 

10  REWIND  2 

NPTS  =  (  1728*N8LOC)/2700 
PRINT  27, NPTS 
READ  TAPE  2,  H 
K  =  1 

TIME  =  -6.5/900. 

M=1 

C  THIS  SECTION  COMPUTES  H  VS  TIME 
DO  60  L=l,900 

H0NE(L)=(H(K)+H(K+3)+H( K+6) )/3. 

HTWC(L )= (H(K+1 )+H(K+4)+H(K+7) )/3. 

LCOUNT=L 

HTHRE( L )=(H{ K+2 )+H(K+5 ) +H( K+8 ) ) /3. 

TIME  =  TIME  +  6.5/900. 

T(L)  =  TIME 
K  =  K  +  3*NPTS 
MAX  =  K+8 

IF(MAX--  1728)  60,60,80 
80  K=1 

IF(M-NBLOC)  85,95,95  ttt  7 

85  READ  TAPE  2,H 
M=M+1 

60  CONTINUE 
95  PRINT  97,LC0UNT 
PRINT  97 -M 

i*t16/100. 

REWIND  2 

-  WRITE  TAPE  2,  HONE 
WRITE  TAPE  2,  HTWO 

_ _ i.ioTTc  TArirr  o  tjTim  c - - - - - -  .  _ _ _ _  -  . _ _ 


r  luis  t:  M  V  -  u  . 

HTW0AV=0. 

HTHREAV=0. 

DO  U8  L=1,900 

HONEAV=HONEA\/  +  HONE  (  L  )  APPENDIX  III 

HTWCAV=HTWOAV+HTWO(L )  - 

48  HTHREAV  =  HTHREAV+HTHRE (L  ) 

H0NEAV=H0NEAV/900. 

HTW0AV=HTW0AV/900. 

HTHREAV=HTHREAV/900. 

PRINT  49tH0NEAVtHTW0AV, HTHREAV 

49  FORMAT! 3E20. 7) 

PRINT  27,NBL0C 

27  FORMAT! 17) 

HONEAV=  IHONEAV  +HTWOAV  +HTHREAV)/3. 

PRINT  26,HONEAV 
26  FORMAT ! E20. 7 ) 

00  90  L=1 ,900 

HONE!L)  =  HONEIL)  -  HONEAV 
HTWC(L)  =  HTWO!L)  -  HONEAV 
90  HTHREIL)  =  HTHREIL)  -  HCNEAV 
DO  1  L=l,12 
1  ITITLEd)  =  8H 

ITITLE! 1 )=8HAN0ERS0N 
ITITLE!2)=8H  BOX  263 
ITITLE!7)=8H  EARTHS 
ITITLE!8)=8HMAGNETIC 
IT1TL£!9)=8H  FIELD  V 
ITITLE! 10)=8HS  TIME 
ITITLE! 1  I )  =  8HT  IN  HRS 
ITITLE ! 12 )=8H  H  GAMMA 
LA8A=4H  HI 
LABB=4H  H2 
LABC=4H  H3 

CALL  DR AW! 900, T, HONE, 1, 0,LABA, ITITLE, 1.0, 10. 0,5, 1,2, 2, 8, 10,0, LA ST) 
CALL  DRAW1900, T,HTW0,2, O.LABB, ITITLE, 1.0, 10. 0,5, 1,2, 2, 8, 10,0, LA ST) 
CALL  DRAW(900,T,HTHRE,3,0,LA8C, ITITLE, 1.0, 10. ,5, 1,2, 2, 8, 10,0, LA ST) 
CALL  SATGRAF (0, 900,T,H0NE,0,0,0,0) 

CALL  SATGRAF (0,900, T,HT WO, 0,0, 0,0) 

CALL  SATGRAF !0, 900, T,HTHRE, 0,0, 0,0) 

C  THIS  SECTION  COMPUTES  SIMPLE  DIFF 
READ  TAPE  2,  HONE 
READ  TAPE  2,  HTWO 
READ  TAPE  2,  HTHRE 
REWIND  2 
DO  1 5  L=1 ,900 
H(L)=HONE(L) 

HONE!L)=H(L)-HTWO(L) 

STCRE=HTWO(L  ) 

HTWO!L)=H(L)-HTHRE(L) 

15  HTHRE(L)=STORE-HTHRE(L) 

DO  16  L=l,12 

16  IT1TLE!L)=8H 
ITITLE! 1 )=8HANDERS0N 
ITITLE!2)=8H  BOX  263 
ITITLE!7)=8HFIRST  DI 
ITITLE!8)=8HFFERENCE 
ITITLE(9)=8H  IN  MAG 
ITITLE! 10)=8HFIELD 
ITITLE! 1 1 )=8HT  IN  HRS 
ITITLE! 12)=8H  H  GAMMA 
LABA=4HDEL1 
LABB=4HDEL2 
LA8C=4HDEL3 

CALL  DRAW! 900, T , HONE, 1,0,LABA, ITITLE, 1.0, 2. 00, 5, 1,2, 2, 8, 10,0, LA  ST) 
CALL  DRAW ! 900, T, HTWO, 2, C,LABR, ITITLE, 1.0, 2. 00, 5, 1,2, 2, 8, 10,0, LA ST) 
CALL  DRAW!900, T, HTHRE, 3,0, L ABC, ITITLE, 1.0, 2. 0,5, 1,2, 2, 8, 10,0, LAST) 
CALL  SATGRAF !0 , 900 , T , HONE , 0 , 0 , 0 , 0 ) 

CALL  SATGRAF ! 0, 900, T, HTWO, 0,0, 0,0) 

CALL  S ATGRAF ! 0, 900, T, HTHRE, 0,0, 0,0) 

70  CONTINUE 
REWIND  3 
STOP 
END 

MACHINE  READ! 1 ARG,2ARG, 3ARG,4ARG) 


*  NOTE  G  =  ADDRESS  OF  ENTRY  TO  TYPEWRITER  SENSE  ROUTINE 

*  P  =  ADDRESS  OF  PROGRAM  CONTROL  ROUTINE  JUMP  ENTRY 

*  THIS  ROUTINE  IS  FOR  PROGRAM  CALL 

L0C!G*31,  P=10) 

*  TYPEWRITER  MESSAGE  CODE  TABLE 


III-3 


I 


I  fc-  I  -r  . 


CON (T12=45040U1 22030220 4B, 

*  114=14112000000000006, 

*  140=45040401301520046, 

»  142  =  041406221416300  IB, 


113=20062204032604266, 
115=45040431 12 1401 20B, 
141=34061401040603016, 
743=20224200000000006) 


»  TAPE  UNIT  ASSIGNMENT  TABLE 

*  ASSIGNMENTS  FOR  CHANNEL  3/4  LIBRARY  TAPE 


CONIKO 
»  K3 
►  K6 


=0  , 
=3320306, 
=5520206, 


K1 

K4 

K7 


=3320106, 

=3320408, 

=5520306, 


K2 

K5 

K8 


=3320206, 
=5520106, 
=5520406 ) 


*  ASSORTED  CONSTANTS 

C0N(R4  =73737373737373736,  R16=0070000000000000B ) 

*  TAPE  READING  ROUTINE 


IX 

SLJ (N) 

SLJ (L+S) 

.  EXIT/ENTRY 

lARG 

ZRO(O) 

ZRO(O) 

.  TAPE  UNIT  NUMBER  ARG. 

2ARG 

ZRO(O) 

ZROIO) 

.  INITIAL  ADDRESS 

3ARG 

ZRO(O) 

ZRO(O) 

,  TERMINAL  ADDRESS 

4ARG 

ZRO(O) 

ZRO(O) 

.  MODE,  1  =  BINARY,  2  =  DC 

LDA7{1 ARG) 

INA(-1 IB) 

• 

SIUl ( 19X) 

AJP3(L+3) 

.  CHECK  FOR  TU  ASSIGNMENT 

ENA(T40) 

SLJ4{  G) 

.  GO  PRINT  ERROR 

SLJ(  P) 

ZRO(O) 

• 

INA( 116) 

STA( KO  ) 

• 

LILl (KO) 

LDA1 (KO) 

• 

ADD7(4ARG) 

SAU  (4X  ) 

.  SET  TU  CONTROL  CODES 

SCL(777B) 

SALI 4X  ) 

.  INITIALIZE 

SAL(5X) 

SAL( 1 1 X) 

• 

INA(3) 

SAU( 8X  ) 

• 

INA(2) 

SAU(9X  ) 

• 

INA( 1 ) 

SAU( 1 IX) 

• 

INA( 1 ) 

SAU( 7X ) 

• 

LRS( 18) 

ENA(O) 

• 

LLS(3) 

SAL ( 2X ) 

• 

SAL(3X) 

SAU ( 12X) 

• 

LLS(39) 

LDQ{ R16) 

• 

STA(KO) 

LDAISX  ) 

• 

SSU(KO) 

STA(5X) 

• 

LDA7 (2ARG) 

SAU( 5X  ) 

LDA(6X) 

SSUIKO  ) 

• 

STA(6X) 

ENI (0) 

• 

2X 

L0A7{3ARG) 

SALIN) 

• 

SAU( 6X) 

SIL2( 1 9X) 

• 

3X 

ENI2{4) 

LILKN) 

• 

4X 

EXF(N) 

EXF7(N) 

.  SELECT  UNIT 

5X 

EXF (N) 

EXF7(N) 

.  ACTIVATE  CHANNEL 

6X 

EXF(N) 

LDA7(2ARG) 

• 

SUB{R4) 

LLS(48) 

• 

INI! (-1 ) 

LOAl (Z ) 

• 

INI  1 ( 1 ) 

SUB(R4) 

• 

7X 

EXF7(N) 

SLJ( 15X) 

.  SENSE  ENDFILE 

8X 

EXF7{N) 

SLJl L+4) 

.  SENSE  PARITY  ERROR 

9X 

EXF7(N) 

SLJ{  10X) 

.  SENSE  LENGTH  ERROR 

3 

QJPl ( 19X) 

AJP( 3X  ) 

.  SENSE  FOR  BAD  RECORD 

SLJ { 19X ) 

ZRO( 0) 

.  EXIT 

IJP2(L+2) 

SLJ( 14X) 

.  SENSE  FIVE  PARITY  ERRORS 

lOX 

IJP2(L+1 ) 

SLJ ( 13X) 

.  SENSE  FIVE  LENGTH  ERRORS 

AJP(3X) 

QJP ( 3X  ) 

.  SENSE  FOR  RAD  RECORD 

llx 

EXF(N) 

EXF7{ N  ) 

.  BACKSPACE 

12X 

LILl (N) 

SLJ(4X) 

.  RETURN  TO  REREAD 

13X 

ENA{T6) 

SLJ (  17X) 

.  SET  TO  INDICATE  TROUBLE 

14X 

ENA(T9) 

SLJ( 17X) 

• 

15X 

ENA(T12) 

SLJ( 18X) 

• 

17X 

SLJ4(  G) 

ZRO( 0) 

.  TROUBLE  EXIT 

AJP( 19X) 

AJP3(3X) 

.  SENSE  ACTION  TO  TAKE 

SLJ( llX) 

ZRO(O) 

• 

18X 

SLJ4(  G) 

ZRO( 0) 

.  TROUBLE  EXIT 

AJP(3X) 

AJP3(3X) 

.  SENSE-  ACTION  TO  TAKE 

SL J (  1 IX) 

ZRO(O) 

• 

19X 

ENI 1 (N) 

ENI2(N ) 

.  RE-STORE  INDEXES 

SL J( IX) 

END 

END 

ZRO(O) 

• 

III-9 

. 5143083E+05 
.5093095E+05 

29 


.5097462E+05 


APPENDIX  r/ 


CIRCaiT  DLiGP.i\?4S 
POLAHIZII'JG  AKTD  DAIiPING  IffiWOP.K 
MAGIffiT'OMETER  PREAMP 
MAGNETOMETER  APTLIFIER 
RADAR  GATE 

CLOCK  GATE  AND  GATE  CONTROL 
STortiAGE  RESET  AND  GATE  DRIVER 
TIMER 

COUNTER  START  AND  RESET 
SIGNAL  GATE  AND  GATE  CONTROL 
CLOCK  FREQUENCY  MULTIPLIER 
COUNTER  FLIP=FLOP  kW)  GATE 
MULTIPLEX  PULSE  GF.IffiRATOR 
I^TIPLSX  BLOCK  DIAGRAM 
AMLOG  LADDER 
OirrPUT  AP1PLIFIERS 

PLUS  AT®  MINUS  12  VOLT  FOVIER  SUPPLY 

PLUS  km  MINUS  h  VOLT  PO'w’ER  SUPPLY 

PLUS  28  VOLT  OSCILUTOR  AND  TBfER  PCV.'ER  SUPPLY 

INTERFACE  DRIVER  AND  AMPLIFIER 

INTERFACE  AND/OR  GATES 

INTERFACE  PXIP»FLOP  AND  IWERTER 
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